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SECTION  I 


INTRODUCTION 


The  design  of  airborne  electro-optical  systems  requires  the  knowledge 

of  angular  vibration  as  well  as  linear  vibration  of  aircraft  structures. 

Just  as  the  linear  vibration  involves  displacements  in  three  coordinate 

axes,  the  angular  vibration  addresses  itself  to  the  rotary  mode  of  internal 

bending  about  each  of  the  three  axes.  In  one  dimension,  for  instance,  if 

linear  vibration  refers  to  the  transverse  deflection,  the  pitching  motion 

expressed  by  the  local  slope  of  linear  (transverse)  displacement  would 

represent  the  corresponding  angular  vibration.  Theoretically,  it  is  possible 

to  formulate  the  angular  vibratory  response  of  aircraft  structures  subject 

to  aerodynamic  and  acoustic  random  excitations.  Although  modem  computer 

structural  analyses  such  as  the  NAGTRAN  may  be  implemented,  they  are  at  present 

limited  to  the  low  frequency  range  of  vibration  of,  particularly,  a complex 

structure.  To  cope  with  the  high  frequency  range  vibration,  it  is  expedient 

to  relate  the  angular  vibration  directly  to  the  linear  vibration  response. 

This  was  the  approach  of  Reference  1.  Using  the  Bernoulli-Euler  beam  as  a 

theoretical  model,  a simple  qualitative  relationship  was  obtained  between  the 

linear  root-mean-squared  (rms)  amplitude  y and  angular  rms  amplitude  0 

rms  rms 

(Equation  (27)  in  Reference  l) 


0 * 1.17  (n/L  ) y , 

rms  J rms  * 


where  L is  the  beam  length.  An  alternate  relation  was  also  derived  in  terms 
of  the  discrete  angular  deflection  which  is  defined  by  the  ratio  of  the 
difference  of  transverse  displacements  y.,(t)  and  y (t)  at  x0  and  x , respectively, 

*-  1 i.  1 


to  the  separation  distance  Ax  ■ x.,  - x ; namely 

<-  A 


0^(0  - fy^U)  - y^t)]/  Ax. 


Here  the  y's  are  measured  by  using  accelerometers.  The  use  of  the  discrete 

or  differential  angle  0^(t)  i»  dictated  by  the  practical  requirement  that  at 

present  there  is  no  inexpensive  off-the-shelf  type  angular  deflection 

transducer  with  sufficiently  high  frequency  response.  By  a parallel 

derivation,  the  relationship  involving  (0.)  and  y which  is  counterpart 

a ms  ms 

to  (1)  becomes 


(0  ) ■ ( /2/  Ax)  y 

A rms  * rms 


It  must  be  pointed  out  that  relations  (l)  and  (3)  are  independent  of  the  local 
coordinates  and  end  conditions.  Note  that  relation  (3)  depends  only  on  the 
separation  distance.  This  came  about  because  0 and  (0.)  and  y are 
the  rms  amplitudes  averaged  over  the  beam.  The  beam-averaging  originally 
introduced  in  Reference  1 is  based  on  the  assumption  that  response  statistics 
should  not  vary  significantly  from  one  point  to  another  on  a complex  structure 
under  the  temporally  and  spatially  homogeneous  random  excitations.  The 
consequence  of  beam-averaging  is  the  loss  of  detailed  spatial  information. 

This  then  enables  us  to  apply  the  linear-angular  rms  amplitude  relations,  (l) 
and  (3),  to  any  structures  as  a qualitative  relation,  provided  L or  Ax  can  be 
specified  depending  on  the  type  of  angular  measurement  technique  employed. 

Based  on  limited  vibration  data  available  at  the  writing  of  Reference  1, 
it  was  concluded  that  by  relations  (l)  and  (3)  the  angular  rms  value  can  be 
predicted  within  ‘20!t  of  the  measurement.  Since  then,  some  more  linear 


2 


and  angular  vibration  data  have  been  obtained  from  the  flight  tests  of  RF-4C 
and  F-15  fighters,  CH-3E  helicopter,  and  B-52  bomber  (References  2-5). 

This  therefore  enables  us  to  re-evaluate  the  prediction  capability  of 
relations  (l)  and  (3)  in  light  of  the  additional  vibration  test  data. 

Moreover,  it  is  now  possible  to  move  a step  further  and  predict  the  spectral 
shape  of  the  angular  power  spectral  density  (psd)  function,  which  is  the 
main  objective  of  this  report. 

Following  the  same  line  of  arguments  as  in  Reference  1,  a general 
relationship  has  been  derived  between  the  linear  and  angular  response  psd's, 
which  states  that  the  linear  psd  falls  off  faster  than  the  angular  by  the 
power  exponent  of  one  (Section  II).  This  linear-angular  psd  relation  was 
first  tested  on  a simply  supported  beam  excited  randomly  at  the  midspan  in 

such  a manner  that  the  acceleration  response  psd  was  constant  over  40  - 3200  Hz. 

-1* 

Consequently,  the  linear  displacement  psd  has  the  power-form  of  f It 

was  then  found  that  the  angular  psd  obeys  the  power-form  of  f~\  as  predicted, 
over  the  same  frequency  range  (Section  III). 

Encouraged  by  the  success  on  test  beam,  an  attempt  was  made  to  predict  the 
spectral  shape  and  magnitude  level  of  the  angular  psd  of  typical  flight  tests 
of  RF-UC  and  F-15  fighters,  CH-3E  helicopter,  and  B-52  bomber  (Section  IV). 

After  splitting  the  linear  vibration  psd  into  the  low  and  high  frequency 
ranges,  the  angular  vibration  psd  has  beer,  predicted  in  two  separate  frequency 
ranges.  First,  the  spectral  shape  of  angular  psd's  can  be  predicted  accurately 
in  the  high  frequency  range  for  all  test  cases,  but  not  always  so  in  the  low 
frequency  range.  Second,  an  accurate  prediction  of  the  magnitude  level  is 
difficult  to  attain  in  that  it  requires  the  correct  spectral  shape  and  rms 
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SECTION  II 

THEORETICAL  ANALYSIS 

We  shall  derive  in  this  section  a simple  relationship  between  the  linear 
and  angular  psd's  which  is  counterpart  to  the  rms  amplitude  relation  (l). 
Although,  as  in  Reference  1,  the  Bernoul  li-Euler  beam  will  be  used  as  a 
theoretical  model,  it  turns  out  that  the  psd  relationship  derived  here  is 
applicable  to  actual  flight  test  data  (Section  IV). 

We  begin  with  the  Bernoul li-Euler  beam  equation 


where  pA  is  the  mass  per  unit  length  of  the  beam,  0 is  the  viscous  damping 
constant,  El  is  the  flexural  rigidity,  and  f(x,t)  is  the  excitation  force 
per  unit  mass  of  the  beam.  To  facilitate  the  analysis,  split  the  forcing 
function  into  the  spatial  part  (q(x))  and  temporal  part  (g(t)); 

f( x,  t)  - q(x)  g(t).  (5) 


a2y 

♦ 2 0 

ay 

EI 

♦ - 

aUy 

= f ( X , t ) , 

(*> 

at2 

at 

oA 

, L 

3x 

i 

The  stochastic  dynamics  of  the  beam  can  be  formulated  based  on  Equations  (1*) 
and  (5)  (References  6 and  7);  the  appropriate  extension  to  angular  vibration 
has  been  carried  out  in  Reference  1. 


1.  Linear  and  Angular  Response  PSD's 

Theoretically,  it  is  convenient  to  consider  the  simply  supported  beam 
of  length  L because  the  undamped  normal  modes  are  sin(nrx/L).  On  the  other 

i 

I 


I 

I 


haiul,  a simply  supported  beam  is  not  as  experimentally  easy  to  simulate  as  ones 

with  fixed  and  free  end  conditions.  In  any  event,  the  simply  supported  end 
condition  is  inconsequential  because  it  turns  out  that  the  predicted  spectral 
shape  of  angular  pad  is  independent  of  a particular  end  condition  within  the 
present  formulation.  The  spatial  behavior  of  y(x,t)  and  q(x)  can  be 
expressed  by  the  normal  modes 

« oo 

y(x,t)  - l y (t)  sin(nnx/L),  q(x)  - l q sin(n*x/L).  (6) 

n-1  n=l  ^ 

'Ihe  solution  of  y^(t)  will  be  obtained  in  terms  of  the  modal  impulse-response 
h (t)  - [(exp(  -at )//j^2  - B^lsinU/)2  - B2},  where 

n n ■*  n 

■ (ns/L)1*  EI/pA.  (7) 

The  transverse  deflection  subject  to  the  zero  initial  displacement  and 
velocity  is 


y(x,t)  * £ sin(nwx/L)  □ / h (t)  g(t-  t)  dT. 

n-1  0 n 


(8) 


For  a small  amplitude  vibration,  the  angular  deflection  may  be  approximated 
by  the  local  derivative;  8(x,t)  - dy/dx. 


0(x,t)  - (w/L)  I n cos(nnx/L)  q_  / h ( t ) g(t  - t)  dT  . (9) 

_ , n - n 

n-1  0 
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Under  the  assumption  that  g(t)  and  q(x)  are  statistically  independent 
and  that  g(t)  is  s stationary  process,  the  variances  of  y and  t)  are 
(liquations  l 10)  and  (12)  in  Keference  l) 


2 

<y“(x)>  » E sin(n"x/L)  sin(knx/L)  <q  q > I , 

n,k-l  ^ k nk 

90 

2 O 

<0  (x)>  ■ (s/L)*'  E nk  cos(n*x/L)  coa(knx/L)  <n  q > 
n,k-l  ^ k 

where  < > denotes  statistical  average.  The  I , reflects  both  the 

nk 

system  characteristics  and  excitation  force 


I,,k  " (l/1*”)  / H ( u>)  H»(u>)  0 (w)  du>  , (11) 

iwv  n k pg 

2 2 -»X 

where  Hn(u>)  ■ ( ♦ i2Bu>  ) is  the  nodal  frequency  response,  and  * 

denotes  the  complex  conjugate.  Further,  0 (u>)  is  the  forcing  psd  related 

SR 

to  the  correlation  K (t)  ■ <g(t  ♦ i)g(t)>  by  the  Fourier  transform  pair; 


K (t)  = ( 1 /Uit ) / G (w)  exp(  iuiT  ) dm, 
gg  _ gg 

CD 

G (w)  ■ 2 / R (t)  exp(-iwT  ) dr  . 

BO  So 


> 


(12) 


The  factor  (M_1  ensures  that  R (0)  - /"  G ( f ) df  with  the  substitution 

gg  0 gg 

of  ii>  ■ 2xf.  That  is,  the  variance  is  the  sum  of  spectral  energies  of  a 
one-sided  psd  (Reference  8). 

When  the  beam  is  assumed  to  be  excited  by  the  spatially  random  q(x) 
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- — j— t-  - :y.~A 


with  a Dirac -delta  spatial  correlation,  Equation  (10)  becomes  simplified 
greatly  because  of 


'W  ■ K V' 


(13) 


where  K is  a constant,  and  6 . • 1 if  n » k and  » 0 if  n 4 k.  It  has  been 

nK 

shown  in  Reference  1 that  Equation  (13)  is  approximately  valid  under  the 
localized  forcing  of  a slightly  damped  structure.  As  a consequence  of 
Equation  (13),  Equation  (10)  reduces  to  a single  sum 


<y2(x)>  * Z sin2(nn  x/L)  K I , 

n-1  nn 


<02(x)>  =»  (*/L)2  Z n2  cos2(nirx/L)  K I. 
n=l 


nn 


(1U) 


Let  us  now  express  the  variances  in  the  spectral  form 


<y2(x)>  = (l/Un)  /"  G (x,w)  du,  < 02 ( x ) > * (l/l*ir)/“  G (x,u>)  du  , ( 15 ) 

— ® yy  — «»  y u 


which  are  defined  analogous  to  Equation  (12).  Inserting  Equation  (15)  and 

I _ s ( 1/U-H ) / | H ( to)  | 2 G (w)  dw  into  Equation  (lU),  one  obtains  by 

nn  — 00  n CfiT 


nn  ' n ' ' • gg 

equating  the  integrands 


G (x,  u»)  * Z sin  (nnx/L)  |H  (u)|2  K G (w), 

yy  n»l  n gg 


(16) 


G (x,  w)  » ( it/h)2  Z n2  cos2(nirx/L)  |H  (w)|2  K G (u>). 

06  n-l  n 
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They  are  the  linear  and  angular  response  psd's  subject  to  the  arbitrary 
but  stationary  forcing  psd. 


2.  Relationship  Between  the  Linear  and  Angular  Response  PSD's 

To  obtain  the  overall  form  of  psd  which  is  independent  of 

the  particular  end  condition  (and  hence  the  mode  shape).  Equation  (l6)  will 

be  averaged  over  the  beam  (Reference  l).  Denoting  the  beam-averaged  psd's 

by  G (u>)  ■ G (x,w)  dx  and  G (u>)  ■ L”1/  L G (x,w)  dx.  Equation  (l6) 

yy  u yy  uu  u uw 

gives  rise  to 


G («)  - l |H  (u>)  I2  (K/2)  G (u>), 
yy  n«l  n gg 

G - { u>)  * (n/L)2  I n2  1 H ( cu ) | 2 (K/2)  G (id), 
60  n-1 


(IT) 


with  the  use  of  L_1/0Lsin2(  nitx/L)dx  * L”1/^  cos2(mrx/L)dx  * 1/2.  By 

taking  the  ratio  of  Equation  (17),  we  find  that  G (w)/G  (u)  is  independent 

yy 

of  the  forcing  psd 


GeeU) 


G (w) 

yy 


(w/L)‘ 


l “ . n2  |H  (u)|2 
n*l  1 n 1 


(18) 


The  analytical  behavior  of  Equation  (18)  can  be  determined  as  follows; 

2 2 
We  have  u>n  3 n from  Equation  (7).  Therefore,  inserting  n = wn/u>^  into 

Equation  (l8)  and  writing  out  the  summations  in  detail  result  in 
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a*eM 


f uijH^fw)  |2  + (J2  |H2(u)|2  + w^lH^tw)!2  ♦ 


= ( tt/l r (i/Wl)x> 


G (u>) 

yy 


( iH1(u)r  + |h2u)t  ♦ iH3(«)r  + ... 


1 


(19) 


Note  that  the  dominant  contribution  of  | Hn( u>)  | 2 * [(u>2  - ui2)2  + U62u>2  j”1 

is  in  a resonance  frequency  band  around  w = u . Therefore,  when  w = id. 

n 1 

the  curly  bracket  { } term  of  Equation  (19)  denoted  by  £ is  approximately 

2 2 

equal  to  w^|H^(u))|  / |H^(w)|  = u>^.  Similarly,  for  w = u>2  we  have  E * 

and  in  general  Z = a>n  for  u * Hence,  the  loci  of  local  maxima 

of  £ increase  linearly  as  w passes  through  the  successive  resonance 
frequencies.  Since  the  density  of  resonance  frequencies  per  octave  frequency 
decade  increases  very  rapidly,  E » u is  a good  approximation  in  the  high 
frequency  range  as  substantiated  by  the  numerical  computation  of  Figure  1. 
Using  this  approximation.  Equation  (19)  in  terms  of  the  frequency  becomes 
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i 

\ 


G08(f)/Vf)  " U /L  V f* 


vhich  is  the  desired  relationship  between  the  linear  and  angular  response  psd's. 
The  ratio  G00(  f ) /Gyy(  f ) consists  of  two  parts,  the  scaling  factor  ( n /L  f^) and 
the  spectral  shape  in  the  power-form  of  f.  Note  that  the  scaling  factor  re- 
flects the  particular  end  condition  and  beam  properties  through  the  primary 
frequency  f^.  On  the  other  hand,  the  spectral  shape  function  is  universal  since 
it  is  completely  independent  of  the  end  condition.  This  is,  perhaps,  not 
surprising  because  the  frequency  relation  (7)  is  asymptotically  valid  for 
large  n under  different  end  conditions  (Reference  9).  Due  to  similarity 
in  the  modal  formulations  of  beam  and  plate,  one  may  expect  the  same  shape 
function  to  be  valid  on  a plate  regardless  of  the  boundary  conditions.  As 
a matter  of  fact,  something  more  is  possible.  It  will  be  shown  in  Section  IV 
that  the  spectral  shape  relation  is  applicable  to  the  high  frequency  part 
of  the  linear  and  euigular  psd's  of  typical  aircraft  flight  test  data. 


Footnote  *)  A similar  attempt  has  been  made  to  predict  the  spectral  shape 
of  Ge0(f)/Gyy(f)  in  Reference  10,  However,  there  is  a noticeable  difference 

in  the  final  results.  Namely,  Equation  (7)  of  Reference  10  which  is  the 
counterpart  of  Equation  (18)  can  be  put  in  the  form 

® oo 

Ge0(u>)/Gyy(w)  * (n/L)  1 1 n H^w)  sin(nw/2)  |^/|  l sin(nn/2)|2.  (t) 

The  derivation  of  Equation  (+)  i3  not  Justified.  Apart  from  the  factor 
sin(nn/2)  which  picks  out  only  the  odd  terms  in  the  summations,  Equation  (t) 
involves  the  sums  of  the  complex  Hn,  whereas  Equation  (18)  sums  up  the 

real  | | , The  overall  difference  shows  up  when  one  compares  Figure  1 

with  Figure  3 of  Reference  10,  which  increases  steadily  with  a slope  of 
about  3/1*  up  to  700  Hz  and  then  falls  off  abruptly  thereafter. 
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3.  Power-Form  Representation  of  the  Angular  PSD 

Let  us  suppose  that  the  linear  pad  has  a simple  power-form  expression 
“it 

Gyy(f)  a Df  f where  D is  a constant  and  k is  the  power  exponent  in  the 
frequency  range  (F.^  ?2) . In  view  of  Equation  (20),  the  angular  psd  would 
then  have  the  form 

G00(f)  - C f'k  * 1,  • (21) 

2 2 

where  C = ( * /L  f^)D.  For  the  aircraft  application,  it  is  difficult  to 
estimate  the  constant  factor  C.  This  is  because  we  do  not  know  a priori 
the  primary  frequency  of  a locally  dominant  mode  at  the  arbitrary  aircraft 
location.  For  this  reason,  we  shall  propose  an  indirect  means  of  determining 
the  constant  factor  by  way  of  the  angular  vibration  amplitude  predicted  by 
Equations  (l)  and  (3).  Deviating  from  the  original  context  (Reference  l), 
let  us  assume  that  relations  (l)  and  (3)  can  be  applied  to  the  rms  values 
computed  over  a frequency  band  as  well  as  the  entire  frequency  range.  Then, 
by  requiring  the  predicted  0^s  over  the  frequency  band  of  (F^,  F^)  be 
identical  to  Jj,  2 G00(f)df,  we  obtain 


C = 


(2  - k) 


rms 


p-k+2  r-k+2 

F1  " *2 


(22) 


Equation  (21)  together  with  constant  (22)  represents  a power-form  representation 
in  ( F^,  F2),  based  on  the  power  exponent  k and  the  angular  0 predicted 
by  the  procedures  of  Reference  1. 
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SECTION  III 


BEAM  EXPERIMENT 

To  verify  Equation  (21)  on  the  very  beam  model  that  it  was  originally 
derived,  a beam  experiment  was  carried  out  to  generate  matching  sets  of 
Linear  and  angular  response  psd's  (Appendix  A).  A simply  supported 
Aluminum  I-beam  of  T2  in  vas  randomly  excited  at  the  midspan  In  such  a manner 
that  a constant  acceleration  response  psd  was  maintained  at  the  midspan 
over  40-  3200  Hz.  Linear  displacement  psd's  were  obtained  by  BBN  501 
accelerometers.  The  discrete  angle  0^(t)  was  obtained  by  differencing  the 
tvo  accelerometer  signals  according  to  formula  (2) 

Having  maintained  a constant  acceleration  psd  at  the  midspan,  the  linear 

_L 

psd's  at  various  locations  on  the  beam  obey  the  power-form  of  f over  the  entire 
frequency  range.  This  is  most  evident  in  G ( f ) taken  almost  at  the  midspan, 

yy 

as  shown  in  Figure  AlOa  in  Appendix  A.  Away  from  the  midspan,  the  linear 

psd's  become  progressively  more  Jagged,  yet  they  can  all  be  represented  by  the 
power-form  of  f“\  We  shall  examine  here  a set  of  the  linear  G^(  f ) and 
differential  angular  G (f)  measured  at  x,  = 25  in  and  x„  = 29  in  on  the 

eA  1 

beam  (and  hence  Ax  = U in).  They  are  given  in  Figures  A8a  and  A8b  in 
Appendix  A,  which  have  been  reproduced  here  in  Figures  2a  and  2b,  respectively. 

A least-squares  power-form  fit  to  the  linear  psd  of  Figure  2a  yields 
G * 5.Cxl0”^f""\  Hence,  the  angular  psd  predicted  by  Equation  (21)  is 

yy 

2 2 — T -3 

G„„  * (n  /L  f,  )f  G =l.l*xl0  f . On  the  other  hand,  by  applying  a least- 

60  1 yy 

squares  straight  line  fit  to  the  angular  psd  of  Figure  2b,  we  obtain 

-7  -3 

Gq6  = 1.0x10  f as  the  measurement.  Although  the  predicted  constant  factor 
is  UoS  higher  than  the  measurement,  what  is  noteworthy  is  the  overall  spectral 
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of  the  angular  psd  that  has  correctly  been  predicted  by  the  power-form 
' over  the  entire  frequency  range. 


SECTION  IV 


AIRCRAFT  FLIGHT  TEST  DATA 

It  has  been  found  in  Section  III  that  relation  (21)  can  predict  the 
spectral  shape  of  angular  psd's  of  the  test  beam  exactly  over  the  entire 
frequency  range,  whereas  the  constant  factor  is  predicted  only  approximately 
within  h0%  of  the  measurement.  It  must  be  noted  that  the  1*0/5  difference 
in  constant  factors  amounts  to  a difference  of  less  than  20J»  in  the  predicted 
and  actual  angular  rms  values,  which  is  well  within  the  prediction  limit  set 
forth  in  Reference  1.  For  the  aircraft  application,  the  constant  factor 
will  be  determined  by  Equation  (22)  instead.  It  is  then  evident  that  an 
accurate  prediction  of  8^  is  prerequisite  for  determining  the  constant 
factor  correctly.  With  this  in  mind,  we  shall  examine  here  the  viability 
of  relation  (21)  as  it  is  applied  to  the  aircraft  flight  test  data  recently 
acquired  by  this  laboratory.  Recall  that  the  linear  psd  (Figure  2a)  of  the 
experimental  beam  was  represented  satisfactorily  by  a single  straight  line 
in  the  log-log  plot  over  the  entire  frequency  range.  In  contrast,  it  has 
been  observed  that  some  of  the  linear  psd's  of  aircraft  flight  tests  (e.g.. 
Figures  5&  and  6a)  can  better  be  represented  by  more  than  one  straight  line. 
Although  there  is  no  theoretical  guideline  except  that  the  frequency  band 
should  not  be  too  small  to  violate  the  high  modal  density  assumption  for 
relation  (20),  it  was  decided  to  represent  the  linear  psd  by  two  straight 
lines  over  all  frequency  ranges  but  the  high  frequency  tail  corresponding 
to  noise.  Admittedly,  the  splitting  of  linear  psd  into  two  parts,  namely, 
the  low  and  high  frequency  ranges,  respectively,  is  quite  arbitrary,  hence 
cannot  be  Justified  on  theoretical  grounds.  In  some  cases  (e.g.,  Figure  5a), 
the  choice  of  the  demarcation  frequency  is  self-evident,  whereas  it  is  not  so 
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in  other  cases  (e.g. , Figure  Ua) . In  any  event,  ve  have  here  decided  to 
represent  all  the  linear  psd's  of  aircraft  flight  tests  by  two  straight  lines 
as  a matter  of  uniformity, but  not  of  necessity.  In  other  words,  the  prediction 
procedures  to  be  outlined  shortly  does  not  depend  on  this  particular 
mode  of  representation.  Perhaps,  the  complete  structural  analysis  of  aircraft 
structures  might  suggest  a demarcation  frequency  based  on  the  modal  density. 
However,  until  such  information  is  available,  the  choice  of  demarcation  frequency 
is  more  or  less  empirical;  hence,  no  further  attempt  will  be  made  to  Justify  it. 
The  prediction  procedures  for  the  angular  psd  are  as  follows; 

(i)  The  linear  psd  is  first  divided  into  the  low  and  high  frequency  parts, 
as  mentioned  before. 

(ii)  The  linear  psd  in  each  frequency  range  is  then  approximated  by  a straight 

• lc 

line  in  the  log- log  plot,G  (f)  = Df  , using  a least-squares  fit.  Of  course, 

yy 

the  straight  lines  in  the  low  and  high  frequency  ranges  do  not  necessarily  meet 
continuously  at  the  demarcation  frequency. 

(iii)  The  linear  yms  is  now  computed  by  integrating  0^(f)  over  the  appropriate 

frequency  range.  To  predict  the  angular  © , either  relation  (l)  or  (3) 

rms 

is  used  depending  on  the  type  of  angular  measurement  employed.  In  other 
words.  Ax  or  L must  be  specified. 

(iv)  Finally,  in  each  frequency  range  the  constant  factor  is  determined  by 
Equation  (22). 

In  this  section,  we  shall  examine  the  typical  flight  test  data  of  RF-UC 
and  F-15  fighters,  CH-3E  helicopter,  and  B-52  bomber.  For  the  RF-Uc,  F-15, 
and  CH-3E  tests,  the  angular  psd  was  obtained  by  the  differential  accelerometer 
technique,  hence  the  prediction  scheme  (3)  is  applicable.  On  the  other  hand. 
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the  angular  differentiating  accelerometer  was  used  to  acquire  the  angular 
psd  of  B-52,  so  that  the  angular  rms  value  will  be  predicted  by  Equation  (l) 
with  a suitably  defined  L,  as  will  be  discussed  later. 

RF-UC  Fighter 

Figure  3a  is  the  vertical  displacement  G (f)  taken  at  the  left  wing  root 
of  a RF-Uc  during  0.9  M cruise;  the  corresponding  pitch  G (f)  with  Ax  = 

Va 

2h  in  is  shown  in  Figure  3b.  (They  are  derived  from  Figures  23  and  150  of 

Reference  2.)  Using  this  set  of  psd's  as  a typical  example,  we  shall 

illustrate  the  procedures  (i)  - (iv).  To  begin  with,  (i)  the  linear  G (f) 

is  divided  into  the  low  frequency  range  of  10  - 50  Hz  and  the  high  frequency 

range  of  50  - 1200  Hz.  As  mentioned  before,  the  spectral  tail  beyond  1200  Hz 

will  be  ignored,  since  it  contains  very  little  spectral  energy.  (ii)  In  each 

frequency  range,  the  linear  psd  is  represented  by  a straight  line  as  shown  in 

l , _3  _l* 

Figure  3a.  The  analytical  expressions  are  G = 6.0x10  f and 

yy 

h -2  -3  5 

Gn  = 1.3x10  f"  * , where  the  superscripts  i and  h refer  to  the  low  and  high 

yy 

frequency  ranges,  respectively.  (iii)  The  linear  rms  values  are 

y = O.l^xlO”  in  and  y = 0. 5^x10”  in.  Hence,  using  Equation  (3)  under 
rms  rms 

t -1* 

Ax  = 2h  in,  the  predicted  angular  rms  values  are  (0.)  = 0.83x10  rad  and 

A rms 

(6A^rms  = rad#  Finally,  the  computed  constant  factors  by 

Equation  (22)  are  = l.UxlO-6  and  Ch  = 5.1*xlO”7. 

Therefore,  the  predicted  angular  psd  is 

Gl  (f)  = l.UxlO"6f‘3,  Gh  (f)  = 5.1*xlO"Tf"2,5.  (23) 

0A0A  6AeA 


IT 


' 

I 


! 

‘‘l 

I I 

[ 


For  comparison,  the  angular  psd  in  Figure  3b  fitted  by  a least-squares 
procedure  yields 


°v4(f) 


l.Uxlo-V1-3.  r,h 


GV  <fl 
A A 


7.0xl0~7f"2*5. 


(2U) 


All  the  relevant  data  are  summarized  in  Table  1,  and  the  predicted  psd 
is  compared  with  the  actual  in  Figure  3b.  It  must  be  noted  that  both 
the  constant  factor  and  power  exponent  have  been  predicted  quite  accurately 
in  the  high  frequency  range.  In  the  low  frequency  range,  however,  the 
predicted  power  exponent  is  considerably  smaller  than  the  actual. 


F-15  Fighter 

Recently,  this  laboratory  has  received  some  preliminary  flight  test  data 
of  the  linear  accelerometers  mounted  on  a bulkhead  (fuselage  station  No.  509) 
near  the  forward  upper  pallet  attachment  points  of  a F-15  (Reference  3). 
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Figure  4a  is  the  lateral  displacement  psd  during  0.9  M flight  under  light 

buffeting;  Figure  4b  is  the  corresponding  roll  psd  derived  from  the  two 

vertically  displaced  accelerometers  under  Ax  * 28.9  in.  Here,  ve  divide 

the  linear  Gyy( f ) into  the  low  frequency  range  of  4 - 20  Hz  and  the  high 

frequency  range  of  20  - 200Hz,  as  shown  in  Figure  4a.  The  relevant  results 

are  summarized  in  Table  I.  In  the  high  frequency  range,  there  is  little 

difference  in  the  predicted  and  measured  as  shown  in  Figure  4b. 

®AeA 

As  in  the  case  of  RF-4C,  prediction  of  the  power  exponent  is  again  low 
for  the  low  frequency  part  of  the  angular  psd. 

CH-3E  Helicopter 

To  extend  the  prediction  scheme  to  yet  another  situation,  we  have  presented 
in  Figure  5 the  linear  and  discrete  angular  pad's  of  a CH-3E  helicopter 
during  40  knots  flight  with  zero  degree  of  roll.  Figure  5a  is  the  lateral 
displacement  psd  taken  at  the  rear  ramp;  the  corresponding  yaw  psd  measured 
by  linear  accelerometers  under  Ax  3 56.5  in  is  shown  in  Figure  5b.  (They 
are  derived  from  Figures  18  and  25  of  Reference  4.)  Unlike  Figures  3 and  4, 
the  linear  and  angular  psd's  of  CH-3E  have  a series  of  spikes  corresponding 
to  the  discrete-frequency  vibration  induced  by  the  main  and  tail  rotor 
blades.  Overlooking  such  spikes,  we  shall  formally  apply  the  prediction 
procedures  by  splitting  the  linear  G ( f ) into  the  low  frequency  range  of 

yy 

4 - 15  Hz  and  the  high  frequency  range  of  15  - 220  Hz.  The  predicted 
angular  psd  is  summarized  in  Table  I and  compared  in  Figure  5b  with  the 
measurement.  The  power  exponents  have  been  predicted  satisfactorily  in 
both  the  low  and  high  frequency  ranges.  On  the  other  hand,  the  predicted 
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magnitude  levels  deviate  from  the  measurement  almost  by  plus  and  minus  6 db 
(based  on  lOxlog [measured  pad  level /predicted  psd  levelj  ) in  the  low  and 
high  frequency  ranges,  respectively.  This  may  be  attributed  to  the  uncertainty 
in  ms  values  when  the  psd's  are  infested  by  spikes  representing  discrete- 
frequency  rotor  resonances. 

B-52  Bomber  j 

Lastly,  ve  shall  examine  a typical  set  of  linear  and  angular  psd's  of 

the  B-52  flight  test  during  initial  take-off.  Figure  6a  is  the  vertical 
displacement  psd  taken  at  the  forward  end  of  an  armor  storage  compartment  near 

the  bulkhead,  fuselage  station  Ho.  1853;  Figure  6b  is  the  corresponding 
angular  psd  of  pitch  component  fluctuations  measured  by  an  angular 
differentiation  accelerometer  (ADA).  (They  are  derived  from  Figures  68  and 
73  of  Reference  5.  ) After  splitting  the  linear  psd  into  the  low  frequency 
range  of  10  - 70  Hz  and  the  high  frequency  range  of  70  - 500  Hz,  the 
prediction  procedures  have  been  repeated  but  with  the  angular  rms  value 
predicted  by  Equation  (l)  under  L =72  in.,  which  is  the  linear  dimension  of  the 
test  platform.  From  the  results  compiled  in  Table  I and  depicted  in  Figure  6b, 
one  finds  that  the  spectral  shape  has  been  predicted  accurately  for  both 
and  G^0.  In  contrast,  the  predicted  constant  factors  are  off  almost 
by  *10  db.  Although  this  cannot  be  explained  away  in  certainty,  it  is 
suspected  that  the  angular  psd  measurement  by  an  ADA  transducer  is  not 
trustworthy.  In  particular,  the  angular  psd  (Figure  6b)  shows  much  higher 
spectral  energy  buildup  in  the  high  frequency  range  than  any  of  the  discrete 
angular  psd's  of  RF-UC,  F-15,  and  CH-3E.  Other  angular  psd's  of  the  B-52 
test  exhibit  a similar  spectral  energy  buildup,  as  has  already  been  observed 
in  another  flight  test  which  also  employed  the  ADA  transducer  (Reference  11). 
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Whether  this  peculiarity  is  simply  a quirk  of  the  ADA  or  a genuine  structural 
behavior  is  yet  to  be  resolved. 

In  summary,  the  angular  psd  of  the  typical  flight  tests  of  RF-UC,  F-15, 
CH-3K,  and  3-52  has  been  predicted  with  a varying  degree  of  success. 

First,  the  spectral  shape  of  angular  psd  can  be  predicted  quite  accurately 
in  the  high  frequency  range,  but  not  always  so  in  the  low  frequency  range. 

This  is  in  a way  expected  because  the  linear-angular  psd  relation  (20)  is 
an  asymptotic  expression  valid  for  the  high  frequency  range.  Second,  it  is 
difficult  to  accurately  predict  the  magnitude  level  of  angular  psd  in  that 
the  correct  power  exponent  and  rms  amplitude  of  the  angular  vibration  must  be 
available.  Last  but  not  least,  it  must  be  pointed  out  that  the  prediction 
procedures  are  subject  to  the  numerical  variations  of  0.1  in  the  power  exponent 
and  10%  in  the  magnitude  level.  This  is  due  to  the  uncertainty  in  estimating 
the  power  exponent  and  constant  factor  of  a straight  line  representation 


of  the  psd  data 


SECTION  V 


OVERALL  ASSESSMENT  OF  THE  PREDICTION  SCHEME 

The  angular  pad  prediction  has  been  carried  out  in  tvo  steps.  In  the 

first  step,  the  exponent  (l-k)  of  power-form  expression  (21)  is  estimated. 

The  power  exponent  has  been  predicted  very  accurately  for  the  high  frequency 

range  of  angular  psd's.  On  the  other  hand,  the  predicted  exponent  deviates 

from  the  actual  as  much  as  1.3  in  the  low  frequency  range.  This  is,  however, 

not  unexpected  because  the  linear-angular  psd  relationship  derived  here 

is  an  asymptotic  expression  valid  for  the  high  frequency. 

The  second  step  involves  estimating  the  constant  factor  by  Equation  (22) 

based  on  the  power  exponent  k and  angular  rms  e . Deviating  from  the 

rms 

original  derivation  in  Reference  1,  it  was  assumed  that  relations  (l)  and 

( 3)  can  be  applied  to  the  rms  values  computed  over  a limited  frequency 

band  as  well  as  the  entire  frequency  range.  This  assumption,  however, 

appears  to  be  a weak  link.  In  the  B-52  case,  for  instance,  the  predicted 

6h  underestimates  the  actual  value  by  threefold,  whereas  the  predicted  0 
rms  rms 

over  the  entire  frequency  range  is  within  50!t  of  the  measurement.  Note, 

however,  that  this  prediction  accuracy  is  worse  than  the  *20/5  limit  claimed 

for  the  total  0 in  Reference  1. 
rms 

Of  course,  the  overall  assessment  of  the  prediction  scheme  must  take 
into  account  both  the  power  exponent  and  constant  factor.  From  Figures  2b, 
3b,  4b,  and  5b,  it  is  concluded  that  in  both  the  low  and  high  frequency 
ranges  the  predicted  angular  psd  lies  within  *10  db  of  the  measurement. 

Though  crude,  the  prediction  scheme  will  be  useful  in  the  preliminary  design 

stage  of  electro-optical  systems  because  it  can  provide  a quick  assessment  of 
the  angular  vibration  environment  prior  to  fabrication.  To  build  up  further 
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confidence,  we  hope  to  test  the  prediction  scheme  on  some  new  flight  test 
data  as  they  soon  become  available  to  us,  perhaps,  with  the  purpose  of 
establishing  certain  confidence  intervals. 

In  retrospect,  the  major  thrust  of  this  report  was  to  propose  a simple 
formula  to  predict  the  angular  vibration  psd  of  aircraft  structures  using 
the  least  amount  of  information,  i.e.,  the  linear  psd  and  a length  scale 
(Ax  or  L).  Consequently,  the  prediction  capability  is  limited.  The 
overall  performance,  however,  is  quite  encouraging  in  that  the  present  work 
j is  Just  an  exploratory  attempt  of  first  order.  Clearly,  the  more  information 

is  utilized,  the  closer  one  should  be  able  to  predict  the  angular  psd. 

Of  course,  the  numerical  structural  analysis  of  an  aircraft  could  yield  the 
detailed  angular  vibration  characterizations  at  every  point  on  the  aircraft 
structure,  subject  to  arbitrary  aerodynamic  and  acoustic  random  excitations. 
Unfortunately,  this  does  not  seem  feasible  at  the  present  state-of-the-art 
of  computer  structural  analysis. 

! 

: 

j 
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Figure  1.  Numerical  Computation  of  E ; 

E = 107  lb/in2,  I = 1.855  in4,  p = 0.131  lb/in 
A = 1.244  in2,  L = 72  in,  6 - w./20 
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Figure  2b.  Angular  Power  Spectral  Density  of  the  Test  Beam 
at  Xj  = 25  in.  and  X?  = 29  in. 

Least-Squares  Approximation:  1.0  x 10’V 

Prediction:  1.4  x 10”7f”3 
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Figure  3a.  Linear  Power  Spectral  Density  of  the 
RF-4C  Flight  Test  (Vertical) 


Least-Squares  Approximation 
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Figure  4a.  Linear  Power  Spectral  Density  of  the 
F-15  Flight  Test  (Lateral) 


Least- Squares  Approximation 
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Figure  4b.  Angular  Power  Spectral  Density  of  the 
F-15  Flight  Test  (Roll) 

Least-Squares  Approximation  ; 

Prediction 
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Figure  5a.  Linear  Power  Spectral  Density  of  the 
CH-3E  Flight  Test  (Lateral) 


Least -Squares  Approximation 
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Figure  5b.  Angular  Power  Spectral  Density 
CH-3E  Flight  Test  (Yaw) 
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Figure  6a.  Linear  Power  Spectral  Density  of  the  B-52 
Flight  Test  (Vertical ) 


Least-Squares  Approximation 
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APPENDIX  A 


BEAM  EXPERIMENT 

This  appendix  describes  the  experiment  of  a simply  supported  Aluminum 
I-beam  to  generate  accurate  linear  and  angular  vibration  data.  The  unique 
feature  of  the  experiment  was  to  excite  the  beam  so  that  a nearly  constant 
acceleration  response  psd  was  produced  on  the  beam.  This  would  facilitate 

the  verification  of  the  linear- angular  psd  relationship  developed  in  the 

. -4 

main  text,  since  the  linear  displacement  psd  will  have  power-form  f 

over  the  entire  frequency  range.  Although  the  initial  intent  was  to 
check  the  accuracy  of  discrete  angular  measurements  by  a prototype 
magnetohvdrodynamic  (MHD)  angular  rate  sensor,  this  has  not  been  success- 
ful due  to  the  sensor  failure.  Therefore,  to  provide  an  independent  check, 
the  beam  experiment  was  numerically  simulated  by  the  NASTRAN  finite-element 
program.  We  shall  begin  by  establishing  a noise  floor  pad  to  assess  the 
accuracy  of  angular  spectral  measurement  by  the  differential  accelerometer 
tecl.  e. 

I.  Noise  Floor  PSD 

At  present,  there  is  no  angular  vibration  transducer  of  the  off-the- 
shelf  type  that  is  as  versatile  and  inexpensively  available  as  the  linear 
accelerometer.  In  practice,  one  therefore  measures  the  angle  by 
differencing  the  signals  of  two  accelerometers  separated  by  a preassigned 
distance.  In  other  words,  approximate  the  true  angle  0(x,tl  by  Equation  (2) 
repeated  here  for  reader's  convenience 

eA(t)  - (y2(t)  - yi(t)]/Ax,  (2) 
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where  v^(t)  and  y.,(t)  are  the  transverse  displacements  at  x^  and  x0, 

respectively,  and  Ax  * x.,  - x^  is  the  separation  distance.  Using  the 

center  location  x * (x,  + x„)/2,  we  have  x,  - x - Ax/2  and 
c 1 2 1 c 

x.,  • x^  + Ax/2;  hence,  the  discrete  or  differential  angle  0^(t)  can  be 

expressed  in  terms  of  x and  Ax. 

c 

Formal  manipulations  of  definition  (2)  will  shed  light  on  the 
limitation  of  discrete  angular  measurements.  First,  obtain  the  variance 
relat ionship 

■-0‘(t'>s  - — — -x  (<yT (t)>  + <y‘(t)'-  - 2<y.  (t)y,(t)N  ).  (23) 

(Axr  1 

The  right-hand  side  is  poorly  conditioned  in  that  it  will  fall  below  the 
2 2 

noise  floor  for  <y^»  3 <jr2>  * <5,15,2>  a8  becomes  smaU.  This 
therefore  imposes  a minimum  Ax  for  the  discrete  angular  measurement. 

N 

To  analyze  this,  let  us  split  v^(t)  into  the  true  signals  v^(t)  and 
the  noises  Oj(t)  (for  i = 1 and  2).  Under  the  usual  assumption  that 
ej(t)  are  statistically  Independent  of  y (t)  and  of  each  other,  one  finds 
that 

<e?(t>>  - <e?(t)'  + — -x-  (<e;(t)>  + <e;(t)>),  (24) 

A A (Ax) 1 1 

where  <0A(t)N  is  given  bv  Equation  (23)  in  which  the  v (t)  are  replaced 
A i 

a 2 2 2 

bv  y.(t).  Further  assuming  <e^(t)>  = <e,(t)>  * <e  (t)>.  Equation  (24) 
i 12 

reduces  to 

<e?(t)>  - <ej(t)>  + — ^“2  < e2(t)-,  (25) 

(Ax) 
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and  Its  psd  analog  Is 


Ge  e <f)  'cH(f,+ 

Va  Va  (Ax) 


2 C (f)  . 
l ee 


(26) 


Because  of  the  factor  (Ax)  " in  Equations  (25)  and  (26),  the  signal-to- 
noise  ratio  would  improve  as  Ax  increases.  Although  this  is  desired. 


Ax  cannot  be  increased  indefinitely  without  deteriorating  rt -solution 

(Reference  Al).  For  better  resolution,  one  must  operate  in  the  opposite 

range  of  small  Ax.  However,  as  Ax  decreases  signal  may  be  overshadowed 

_ 2 

by  the  noise  term  which  increases  by  the  factor  (Ax)  . This  therefore 
sets  a minimum  Ax  for  the  acceptable  signal-to-noise. 


Now,  to  establish  a reference  psd  Gg  g (f)  corresponding  to  a prescribed 

A A 9 

signal-to-noise  ratio  (S/N) , let  the  ratio  G.,  Q (f)/(2G  (f)/(Ax)  ) be 

2 °A0A  ee 

(S/N)4"  (and  not  (S/N)  because  the  psd  refers  to  the  spectral  energy). 

Then,  by  rearrangement 


(f) 


2(S/N)2 

(Ax)2 


C (f>. 

ee 


(27) 


Typically,  the  noise  floor  of  accelerometers  has  a flat  psd  in  the  noise 

-4 

acceleration;  hence  the  noise  psd  is  G (f)  = c,  f , where  c.  is  a constant 

ee  1 1 

to  be  determined  later.  In  view  of  this.  Equation  (27)  becomes 


2 c^ (S/N) ‘ 


G0  0 - 2 
a A (Ax) Z 


r~4 


(28) 


For  a given  (S/N),  Equation  (28)  therefore  defines  a reference  psd;  angular 
spectral  measurements  falling  below  it  should  be  rejected. 
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II.  Experimental  Beam  Test 


1.  The  Test  Set-Up 

A simply  supported  Aluminum  I-beam  was  randomly  excited  at  the  mid- 
span in  such  a manner  that  a constant  acceleration  response  psd  over 
40-3200  Hz  was  maintained  at  the  midspan.  For  the  simply  supported  end 
condition,  the  beam  was  fastened  to  the  supporting  bases  by  a pin  support 
mechanism  consisting  of  precision  needle  bearings  and  high-temperature 
hardened  steel  shafts.  The  supporting  bases  consisting  of  a 12"  high  and 
8"  diameter  solid  Aluminum  pedestal,  were  bolted  to  a massive  vibration 
test  table  as  shown  in  Figure  Al.  Although  the  rotation  centers  of 
bearings  were  not  on  the  central  bending  axis  of  the  beam,  this  would  not 
result  in  significant  inplane  stresses  because  the  transverse  displacement 
was  kept  small  during  the  experiment.  The  effective  beam  length  between 
the  bearing  centers  was  72  in. 

A 25  lb.  electrodvnamic  shaker  (Ling  Model  60  was  attached  to  the 
bottom  flange  at  the  midspan  of  the  beam  via  a fusible  link,  thereby 
transmitting  only  the  vertical  force  and  no  moments  to  the  beam.  To 
maintain  t constant  acceleration  response  psd  at  the  midspan,  the  shaker 
power  amplifier  was  monitored  and  controlled  by  a Time/Data  Vibration 
Control  and  Analysis  computer.  This  control  computer  compares  the  actual 
acceleration  psd  at  the  midspan  with  a prescribed  one  and  then  minimizes 
the  difference.  Figure  A2  shows  a typical  acceleration  psd  that  the  control 
computer  has  maintained  constant  over  40  - 3200  Hz.  The  dotted  lines 
represent  a ‘6  db  band  about  the  control  reference  psd. 
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Test  Procedures 


Linear  displacements  were  measured  by  BUN  501  accelerometers.  The 
differential  angle  was  obtained  by  electronically  differencing  the 
signals  of  two  accelerometers  and  then  dividing  by  the  separation 
distance  as  indicated  by  Kquation  (2).  Accelerometer  signals  were 
amplified  by  Intech  amplifiers  and  recorded  on  I'M  tapes  (15  ips,  27K  Hz 
center  frequency,  and  5K  Hz  bandwidth).  Only  the  even  channels  of  a 
recorder  (Honeywell  9600)  were  used  to  minimize  the  possible  phase  error 
introduced  by  the  tape  heads.  Both  the  Individual  accelerometer  and 
electronically  differenced  signals  were  recorded.  Prior  to  data  acquisition, 
tlte  two  accelerometers  used  for  discrete  angular  measurements  were  balanced 
to  ensure  the  phase  and  magnitude  compatibility.  This  was  done  bv 
applying  Ig  excitation  at  60  Hz  to  both  accelerometers  placed  on  the  lg 
calibration  table  (B  & K TV 6 106) . After  having  amplified  each  accelerometer 
signal,  their  difference  was  displayed  on  an  oscilloscope.  A balance  pot 
placed  in  line  with  one  of  the  accelerometers  before  amplification  was  then 
adjusted  until  a null  difference  was  obtained.  in  addition,  the  signal 
conditioning  electronics  was  calibrated  by  recording  a 1000  Hz  sine  wave 
at  40  db  amplification  on  each  channel.  This  amplification  level  was  t ho 
one  used  in  the  experiment.  Figure  A3  is  the  schematic  diagram  of 
instrumentation  set-up. 

The  test  procedures  may  be  summarized: 

(1)  A prescribed  acceleration  psd  over  40-3200  Hz  was  entered  into  the 
control  computer  system. 
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(2)  For  the  discrete  angular  measurement,  two  accelerometers  were  placed 

on  the  beam  separated  by  a predetermined  Ax  at  the  location  x . A MHD 

angular  rate  sensor  was  mounted  at  x * 27".  Although  it  failed  to 

c 

operate  properly,  it  was  kept  in  its  place  during  the  experiment. 

(3)  The  shaker  feedback  control  system  was  activated.  After  a steady- 
state  excitation  had  been  attained,  the  accelerometer  signals  including  the 
differenced  signal  were  simultaneously  recorded. 

3.  Summary  of  Test  Results 

All  data  were  recovered  through  the  FM  playback  equipment,  filtered 

for  anti-aliasing,  and  digitized  onto  digital  magnetic  tape  via  an  1T1 

Analog- to-Digi tal  converter.  Psd's  were  computed  with  a 9.8  Hz  bandwidth 

using  the  fast  Fourier  transform  on  a Raytheon  704  computer.  The  accelera- 

-4 

tion  psd's  are  multiplied  by  the  factor  ( *?tt  f)  to  yield  the  corresponding 
displacement  psd,  which  is  permissible  for  a stationary  random  process. 

For  the  first  five  tests.  Ax  was  varied  from  4 in  to  12  in  with  a 2 in 

increment  at  x - 27  in.  Four  additional  tests  were  then  performed  with 
c 

Ax  ■ 4 in  but  by  varying  xc«  They  are  designated  by  the  test  numbers  in 


Table  A1 


Test  No 

X 

c 

(in) 

Ax 

| (in) 

Linear  G (f) 

ylYl 

at  x,  * x - Ax/2 

1 c 

Angular  G„  g (f) 

°A  A 

1 

19 

4 

Figure  A7a 

Figure  A7b 

2 

27 

4 

Figure  A8a 

Figure  A8b 

3 

27 

6 

Figure  A8e 

4 

27 

8 

Figure  A8d 

5 

27 

10 

Figure  A8e 

6 

27 

12 

Figure  A8f 

7 

32 

4 

Figure  A9a 

Figure  A9b 

8 

36 

4 

Figure  AlOa 

Figure  AlOb 

9 

50 

4 

1 

Figure  Alla 

Figure  Allb 

Table  Al.  Summary  of  Test  Results 


The  accelerometer  locations  defined  in  terms  of  the  center  location  x 

c 

and  separation  distance  Ax  are  shown  in  Figure  A4. 


Included  in  this  report  are  only  some  of  the  linear  G (f)  at 

V V 

ylyl 

x^,  since  they  are  representative  of  the  others.  For  the  same  reason, 

none  of  the  linear  G (f)  at  x„(=  x + Ax/2)  are  explicitly  shown  here. 

y2y2  2 C 

Note  that  x of  Test  No.  1 is  the  farthest  from  the  midspan.  The  x 's 

c c 

of  Test  Nos.  2,  7.  and  8 get  closer  to  the  midspan,  until  x is  the 

c 

midspan  for  Test  No.  8.  The  test  results  are  first  compared  with  the 
NASTRAN  prediction  in  the  next  section  and  then  discussed  further  in 
Section  IV  of  this  appendix. 
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III.  Finite-Element  Simulation  of  the  Beam  Experiment 
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As  mentioned  before,  it  was  originally  intended  to  check  the 
accuracy  of  discrete  angular  measurements  by  a MHD  angular  rate  sensor. 
However,  this  prototype  sensor  had  failed  during  the  experiment,  hence 
it  could  not  be  used  for  comparison.  To  provide  an  independent  check, 
it  was  therefore  decided  to  simulate  the  beam  experiment  by  a finite- 
element  numerical  method. 

1.  The  NASTRAN  Model 

The  Level  15  NASTRAN  finite-element  code  was  used  to  numerically 
simulate  the  experimental  beam  (References  A2  and  A3).  The  finite-element 
model  whicti  includes  the  attached  masses  of  shaker  armature  and  MHD  sensor, 
consists  of  15  NASTRAN  CBAR  elements,  as  shown  in  Figure  A4.  There  are  in 
all  30  degrees  of  freedom;  14  for  the  vertical  modes  and  16  for  the 
rotational  modes  along  the  axes  normal  to  the  y-  and  x-coordinates . The 
end  nodes  were  vertically  constrained  to  represent  the  simply  supported 
ends.  The  remaining  node  points  correspond  either  to  the  accelerometer 
locations  or  the  center  locations  for  discrete  angle.  The  distance  between 
any  adjacent  nodes  is  less  than  6",  thus  permitting  excitation  of  at  least 
the  first  six  bending  modes,  so  as  to  correctly  represent  the  frequency 
range  up  to  3200  Hz.  The  shaker  armature  (0.79  lb.)  and  the  MHD  sensor 
with  supporting  brackets  (1.18  lb.)  were  incorporated  as  a point  mass  into 
nodes  No.  9 and  No.  7,  respectively.  The  appropriate  material  constants 
such  as  E = 10^  Ib/in",  I = 1.855  ln\  p = 0.131  lb/in\  and  A = 1.244  in", 
and  the  structural  damping  constants  (determined  experimentally  in 
Appendix  B)  were  introduced  into  the  NASTRAN  program,  a brief  outline  of 
which  can  be  found  in  Appendix  C. 
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2.  The  Forcing  PSD 


To  simulate  the  experiment,  it  is  necessary  to  reconstruct  an 
equivalent  forcing  psd  at  the  midspan  that  would  have  produced  the 
constant  acceleration  response  psd.  Reference  A4  describes  how  this  can 
be  done.  First,  the  finite-element  model  is  loaded  at  the  midspan  with  a 
constant  lg  force  over  40-3200Hz,  and  then  the  response  acceleration  psd 
is  determined  at  the  midspan.  Clearly,  the  acceleration  psd  will  not  be 
flat,  nor  will  it  be  of  the  same  magnitude  as  the  experiment.  Therefore, 
the  magnitude  of  the  forcing  psd  in  each  frequency  band  must  be  adjusted  so 
that  the  computed  response  acceleration  of  the  NASTRAN  model  will  match 
the  experiment.  Figure  A5  is  a typical  forcing  psd  necessary  to  produce 
a constant  acceleration  response  psd  shown  in  Figure  AO,  in  which  the 
NASTRAN  and  experimental  acceleration  psd's  are  compared  with  good  agree- 
ment . 

3.  Comparison  of  the  NASTRAN  results  with  Experimental  Data 

It  is  observed  from  Figures  A7-A11  that  NASTRAN  can  qualitatively 
predict  the  beam  experiment,  although  some  quantitative  differences  exist 
in  the  resonance  frequency  and  its  peak  height.  In  general,  the  agree- 
ment between  the  NASTRAN  result  and  experiment  is  much  better  for  G _ (f) 

ylVl 

than  for  (f).  This  is  not  at  all  unexpected  because  differential 

o.o. 

A A 

angular  measurements  are  vulnerable  to  the  signal  differencing  error. 
Clearly,  the  worst  comparison  is  observed  in  Figure  A10  for  the  angular  psd 
at  the  midspan;  the  NASTRAN  C^  ^ (f)  is  much  lower  than  the  experimental 
In  the  very  low  frequency  range.  This  indicates  that  the  experimental 
beam  deviates  from  the  simply  supported  end  condition  and  midspan  local 
excitation.  According  to  the  theoretical  analysis  (Equation  16),  the 
angular  psd  is  zero  at  the  midspan  because  only  the  odd  modes  can  survive 

ue 


under  the  local  midspan  forcing.  In  reality,  however,  the  spatial 
function  is  not  exactly  cos~’(mix/L)  due  to  the  added  masses  of  shaker 
armature  and  MUD  sensor,  which  would  certainly  result  in  unsymmetric  mode 

shapes.  Nevertheless,  the  low  level  of  C (f)  predicted  hv  NASTRAN 

A A 

appears  to  he  consistent  with  the  unsymmetric  mode  considerations. 


VI.  Discussion  of  the  Experimental  Results 


The  discrete  angular  signal  was  generated  by  instantaneously 
differencing  the  signals  of  two  accelerometers  separated  by  Ax. 

To  detect  balancing  errors,  the  two  accelerometer  signals,  replayed  from 


FM  tape,  were  digitally  differenced  to  compute  the  differential  angle. 
As  anticipated,  there  is  no  difference  in  the  two  angular  psd  results. 


This  is  seen  by  comparing  the  angular  psd  of  Figure  A8b  with  Figure  A12 
which  is  the  angular  psd  obtained  by  digitally  differencing  the  recorded 
accelerometer  signals.  In  this  section,  we  shall  first  assess  the  accuracy 
of  discrete  angular  measurements  in  terms  of  the  noise  floor  psd  and  then 
examine  the  spatial  variation  of  the  linear  and  angular  rms  amplitudes. 


1.  Noise  Floor  PSD 

The  angular  psd  will  be  compared  with  the  reference  noise  psd  given 
by  Equation  (28).  The  constant  c^  in  that  equation  can  be  evaluated  as 
follows:  The  noise  acceleration  psd  of  BBN  501  accelerometer  is  flat 

_ 6 2 

over  5-25K  Hz,  with  the  total  spectral  energy  of  2.25x10  g . From  the 

“6  2 

ordinate  of  the  noise  acceleration  psd,  2.25x10  g /25.000  Hz  or 

1.34x10  ^ (in/sec^) “/Hz,  one  finds  that  c^  = 1.34x10  ^/(2tt)^  = 8.6x10  ^ 

2 

in  / Hz.  Inserting  this  into  Equation  (28)  yields 


Wf) 

A A 


1.72x10  8 (S/N)2  -4 

2 1 

(Axr 


For  (S/N)  = 10  and  10“"  corresponding  respectively  to 

10  db  and  20  db,  the  reference  psd's  are  drawn  in  all  figures  for  the 

angular  psd  (Figures  A7b,  A8b-A8f,  A9b,  AlOb,  and  Allb).  Since  the 


accelerometer  noise  is  but  one  of  many  noise  sources,  we  shall  consider 


the  (f)  for  (S/N)  - 10  db  as  a total  noise  floor.  Note  that  all 

c A A 

angular  pad's  lie  above  the  total  noise  floor  in  all  frequency  ranges; 

hence,  the  angular  spectral  measurements  are  accepted  with  sufficient 

confidence.  An  angular  psd  can,  however,  exceed  the  G.  . (f)  for 

lA  A 

(8/N)  ■ 20  dh  onlv  In  the  high  frequency  range. 

For  the  effect  of  variable  Ax,  let  us  examine  the  angular  pad's  in 
Figures  A8h-A8f.  It  Is  seen  that  the  angular  psd  moves  up  higher  above 

the  reference  U.  . (f)  as  Ax  increases.  This  means  that  we  are  more 

1 A A 

confident  in  the  angular  psd  measurement  as  Ax  becomes  large.  In 
particular,  the  angular  psd  for  Ax  =12  in (Figure  A8f)  lies  well  above  the 

G.  t ft)  for  (S/N)  “ 20  db.  One  cannot,  however,  increase  Ax  indefinitely. 

A A 

I'liis  is  because  for  some  higher-order  normal  modes  there  will  occur  more 
than  one  node  between  the  two  accelerometers,  hence  definition  (2) 
becomes  meaningless  for  such  normal  modes.  A way  of  relating  a maximum 
Ax  with  the  upper  frequency  limit  f*  has  been  suggested  in  Reference  Al, 


As  seen  from  Figures  A8b-A8f  the  upper  frequency  limit  decreases  from 

9146  11*  to  1016  Hz  as  Ax  increases  from  4 in  to  12  in. 

2.  The  1. inear  and  Angular  Vibration  RMS  Amplitudes 

For  the  overall  measure  of  vibration  energy,  the  psd's  are  integrated 
upwards  from  the  frequency  of  40  Hi*,  to  compute  the  rms  values  which  are 
summarized  In  Table  A2.  The  low  frequency  limit  is  well  below  the  primary 
frequency  and  high  enough  to  block  out  the  background  noise.  In  Columns  1 
and  11,  the  rms  values  of  v(Xj)  and  v(x.,)  are  given  for  Test  Nos.  1-8. 


1*9 
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The  arithmetic  means  of  Columns  I and  II  are  shown  in  Column  III,  and 
the  linear  vibration  rms  amplitudes  predicted  by  NASTRAN  are  given  in 
Column  IV.  The  values  of  Columns  III  and  IV  are  plotted  in  Figure  A13 


as  a function  of  x.  Just  as  sin~(n  x/L)  of  the  theoretical  beam 
(Equation  14),  one  finds  that  the  NASTRAN  model  predicts  monotonically 
increasing  linear  rms  values  towards  the  midspan.  Except  the  irratically 
high  value  at  x =27  in, the  experimental  rms  values  also  seem  to  increase 
towards  the  midspan. 


Next,  the  experimental  angular  rms  values  are  given  in  Column  V,  and 

the  computed  values  by  NASTRAN  at  the  center  location  are  shown  in  Column  VI , 

Figure  A13  compares  the  experimental  (0.)  for  Ax  = 4 in  with  the  NASTRAN 

A rms 

results  to  show  spatial  distribution  on  the  beam.  Although  (®^)rms  falls 
off  towards  the  midspan,  it  does  not  become  zero  as  theoretically  predicted 
by  Equation  14  under  the  localized  midspan  excitation.  It  has  been  shown 
that  the  beam-averaged  (0.)  is  related  to  the  beam-averaged  v bv 
(Equation  32  in  Reference  1 of  the  main  text) 


(0.) 

A rms 


r 2 sin(irAx/2L) , 

^ Ax  J^rms' 


Using  the  average  of  Column  III,  yrmg  = 0.19x10  , Equation  (30)  for  Ax  - 

4 in  gives  (0.)  = 0.83x10  ^ rad.  which  is  larger  than  the  local 

° A rms  ° 


measurements  as  shown  in  Figure  A13. 


I 


Table  A2.  The  Linear  and  Angular  Vibration  RMS  Values. 


V.  Conclusions 


To  verify  the  relationship  between  the  linear  and  angular  response 
psd's,  a simply  supported  beam  was  excited  at  the  midspan  to  generate 
a constant  psd  of  linear  acceleration  response  on  the  beam.  The  overall 
consistency  of  linear  and  angular  spectral  measurements  was  demonstrated 
by  numerically  simulating  the  beam  experiment  with  a NASTRAN  finite- 
element  model  which  has  reproduced  the  spectral  details  with  sufficient 
realism. 

Some  pertinent  conclusions  may  be  summarized  as  follows: 

(1)  Angular  spectral  measurements  are  not  affected  by  whether  the 
accelerometer  signals  are  instantaneously  differenced  at  the  time  of  data 
acquisition  or  digitally  differenced  later  from  the  recorded  signals. 

This,  of  course,  presupposes  the  maintenance  of  amplitude  and  phase 
calibrations  of  the  accelerometers. 

(2)  The  angular  spectral  measurements  are  well  above  a noise  floor  psd 
which  is  defined  as  ten  times  the  noise  floor  of  a typical  accelerometer. 
Hence,  the  differential  angular  psd's  are  accepted  with  confidence. 

(3)  From  the  local  rms  values  of  linear  and  angular  vibration,  it  is 
inferred  that  the  simply  supported  end  condition  and  the  localized  midspan 
excitation  are  not  observed  by  the  experimental  beam.  Certainly,  the 
added  masses  of  shaker  armature  and  MHD  sensor  further  contribute  to  the 
divergence  of  experiment  from  the  theoretical  beam  analyzed  in  Section  II 
of  the  main  text. 
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MHD  Sensor 


Accelerometer  #1  i Accelerometer  #2 


Figure  A3.  Schematic  of  Instrumentation  Setup 
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A4.  Accelerometer  and  NASTRAN  Model  Node 
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Figure  A6. 


Response  Acceleration  Power  Spectral  Density 
at  the  Midspan  of  the  Beam 

Measured  ; 

NASTRAN 
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Figure  A7a.  Linear  Power  Spectral  Density  of  Test  No.  1 

Measured  ; 

NASTRAN 
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Figure  A8a.  Linear  Power  Spectral  Density  of  Test  No.  2 


Measured 


NASTRAN 
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FREQUENCY  - HZ  DELTA  F = 9.7656 

Flour?  ARc.  Angular  Power  Spectral  Density  of 
Test  No.  3,  AX  = 6 in.,  Xc  = 27  in. 


— Measured  ; 
. . NASTRAN 
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FREQUENCY  - HZ  pelts  F = 9.7656 


Figure  A8e.  Angular  Power  Spectral  Density  of 

Test  No.  5,  AX  = 10  in. , Xc  = 27  in. 

Measured  ; 

NASTRAN 
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Figure  A8f.  Angular  Power  Spectral  Density  of 

Test  No.  6,  AX  = 12  in.,  Xc  = 27  in. 


Measured  i 
NASTRAN 
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FREQUENCY  - HZ  DELTA  F = 9.7656 


Figure  AlOb.  Angular  Power  Spectral  Density  of 
Test  No.  8,  AX  = 4 in.,  Xc  = 36  in. 

Measured  ; 


NASTRAN 


Figure  AlOb.  Angular  Power  Spectral  Density  of 
Test  No.  8,  AX  = 4 in.,  XQ  = 36  in. 

Measured  ; 


NASTRAN 
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FREQUENCY  - HZ  DELTA  F = 9. 

Figure  Alla.  Linear  Power  Spectral  Density  of  Test  No.  9 

Measured  ; 

NASTRAN 


ZM/ros  awyjQcd 


Figure  Allb.  Angular  Power  Spectral  Density  of 
Test  No.  9,  AX  = 4 in.,  Xc  = 50  in. 

Measured  ; 


NASTRAN 
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Figure  A12.  Digitally  Differenced  Angular  Power  Spectral 
Density  of  Test  No.  2 
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ANGULAR  RMS  AMPLITUDE 
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APPENDIX  B 

THE  DETERMINATION  OF  STRUCTURAL  DAMPING  CONSTANTS 

As  an  input  to  the  NASTRAN  code,  the  structural  damping  constants 
of  the  simply  supported  beam  were  determined  by  a Time/Data  vibration 
Control  and  Analysis  computer  via  the  Modal  Analysis  and  Modeling  Systems 
(MODAMS)  software  package.  The  MODAMS  program  extracts  viscous  damping 
information  from  the  transfer  function  observed  at  a point  on  the  structure 
sibject  to  the  impact  excitation  at  another  point.  The  equivalent  viscous 
damping  coefficients  estimated  by  MODAMS  for  the  first  five  resonance 
frequencies  are  listed  in  Table  Bl: 


Mode  No. 

Resonance 

frequency 

Viscous  damping 
coefficient 

1 

76.9 

0.0374 

2 

265 

0.014 

3 

615 

— 

4 

1302 

0.0215 

5 

1421 

0.00055 

Table  Bl.  Viscous  damping  coefficients 


The  MODAMS  program  failed  to  produce  the  damping  value  of  the  third  mode 
and  predicted  the  fourth  resonance  frequency,  which  is  too  large.  A typical 
output  of  MODAMS  is  shown  in  Figure  Bl, which  depicts  a peak  of  the  transfer 
function  magnitude  about  the  resonance  frequency  of  265  Hz. 
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The  algorithm  of  MODAMS  generates  several  damping  values  from  a least- 
squares  fit  of  the  tranfer  function  magnitude  as  shown  in  the  figure. 

One  then  chooses  the  frequency  (265  Hz)  and  damping  (0.014)  values 
corresponding  to  the  maximum  amplitude  (49.48).  The  remaining  pairs  of 
frequency  and  damping  in  Table  B1  were  obtained  in  a similar  manner. 
Interested  readers  are  referred  to  the  user's  manuals  for  Time/Data 
computer  (1050  Fast  Meadow  Circle,  Palo  Alto,  California)  and  for 
MODAMS  software  (Structural  Dynamics  Research  Corp.,  Cincinnati,  Ohio). 
Further  information  may  be  found  in  "Modal  Survey  Activity  via  Frequency 
Response  Functions,"  by  A.  Klosterman  and  R.  Zimmerman,  SAE  National 
Aerospace  Engineering  and  Manufacturing  Meeting,  November  1975. 

Twice  the  viscous  damping  values  of  Table  B1  were  introduced  into 
the  NASTRAN  program  as  the  structural  damping  coefficients  as  a function 
of  resonance  frequency.  The  input  card  TABDMP1  also  calls  for  the 
structural  damping  values  of  the  lowest  (40  Hz)  and  highest  (3200  Hz) 
frequency  limits.  In  the  absence  of  better  information,  we  have  simply 
assigned  2 x 0.0374  and  2 x 0.00055  as  the  structural  damping  at  40  Hz 
and  3200  Hz,  respectively. 


78 


I 


APPENDIX  C 

THE  NASTRAN  FINITE-ELEMENT  COMPUTER  PROGRAM 

The  NASTRAN  finite-element  code  was  used  because  of  its  extensive 
application  in  the  dynamic  response  analyses  of  aircraft  structures. 

Of  the  many  NASTRAN  solution  formats,  we  have  used  the  NASTRAN  solution 
11*0  for  the  present  stationary  random  dynamic  analysis.  The  NASTRAN 
computer  listings  of  the  test  beam  are  shown  in  Tables  Cl  - C4. 

Tables  Cl  and  C2  are  respectively  the  executive  and  input/output  controls, 
and  Tables  C3  and  C4  are  the  finite-element  model  parameters  and 
prescribed  forcing  psd  input  at  the  midspan  of  the  beam,  respectively. 
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TABLE  C2.  INPUT-OUTPUT  CONTROL 
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To  test  the  accuracy  of  a NASTRAN  model  for  the  experimental  beam, 
the  first  six  resonance  frequencies  are  compared  with  the  experimental 
values  in  Table  C5, which  shows  good  agreement  between  them. 


Mode  No . 

NASTRAN 

Experimental 

1 

55.52 

58.54 

2 

243.36 

253.90 

3 

536.95 

449.21 

4 

939.50 

1044.92 

5 

1478.70 

1503.90 

6 

2207.60 

2109.37 

Table  C5.  Comparison  of  the  NASTRAN  and 

experimental  resonance  frequencies. 
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